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ABSTRACT 


The  central  results  concern  the  initial-value  problem  for  ut  *  + 

u(1-u)(u-a)  ,  -L  <  x  <  L  and  t  >  0,  under  Dirichlet  conditions  at  x  =  ±  L. 
Here  m  >  1  and  0  <  a  <  (m+1  )/(m+3)  .  The  equilibrium  solutions  of  this 
problem  are  determined  for  each  L  >  0  and  it  is  shown  that  the  ur limit  set 
atfUg)  of  an  initial  datum  uQ  with  values  in  [0,1]  consists  of  a 
connected  set  of  equilibria.  This  is  used  to  determine  some  domains  of 
attraction  of  isolated  equilibria.  A  novel  feature  of  the  results  is  that  for 
large  L  there  are  multiple  parameter  families  of  equilibria. 


A  second  part  of  the  paper  gives  a  self-contained  development  of 
existence,  uniqueness,  maximum  principles,  and  continuous  dependence  on  data 
for  more  general  equations  ut  =  n<u)xx  +  f(  u)  •  The  results  are  employed  in 
proofs  of  some  of  the  theorems  referred  to  above. 


Interest  in  these  questions  is  stimulated  by  the  occurrence  of  such 
models  in  science,  e.g.  in  fluid  flow  in  porous  media  and  biology. 


AMS  (MOS)  Subject  Classifications:  35K55,  35K65 

Key  Words s  nonlinear  diffusion,  biological  models,  asymptotic  behaviour, 
equilibrium  solutions,  flow  in  porous  media. 
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STABILIZATION  OF  SOLUTIONS  OF  A  DEGENERATE  NONLINEAR 


DIFFUSION  PROBLEM 

Donald  Aronson,  Michael  G.  Crandall  and  L.  A.  Peletier 

Introduction 

In  this  paper  we  are  primarily  concerned  with  the  large  time  behaviour  of  nonnegative 
solutions  of  the  initial-boundary  value  problem 

fut  -  (um)xx  +  in  x  R+  # 

u(±L,t)  -  0  in  R*  , 

^  u(x,0 )  »  uQ(x)  in  , 

where  m  >  1  is  a  parameter/  f  is  locally  Lipschitr  continuous,  f(0)  *  0,  and  uQ  i3 
bounded.  Problems  of  this  form  arise  in  a  number  of  areas  of  science j  for  instance,  in 
models  for  gas  or  fluid  flow  in  porous  media  [2]  and  for  the  spread  of  certain  biological 
populations  [13,  16] • 

This  paper  is  divided  into  two  part,  in  Part  I  we  consider  what  may  be  called  the 
motivating  example.  Problem  I*,  which  consists  of  Problem  I  with  the  special  choice 
(1)  f(u)  -  u(1-u)(u-a) 

for  suitably  restricted  parameters  a.  We  begin  by  describing  in  detail  the  set  E  =  E(L) 
of  nonnegative  equilibrium  solutions  of  Problem  I*.  Clearly  E(L)  contains  the  trivial 
solution  u  *  0  for  all  L  >  0.  Write 

E*(L)  -  (L) \  {0 }  . 

In  the  description  of  E#(L)  there  are  two  critical  parameter  values  Lg  and  Lj  with 
0  <  Lq  <  L^  <  We  show  that: 
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(i)  E*(L)  =■  $  for  0  <  L  <  L0 
(il)  E  (Lq)  consists  of  one  isolated  positive  solution 
(iii)  For  Lq  <  L  <  L^,  E  (L)  consists  of  two  isolated  positive  solutions,  p  and  q, 
with  p  <  q  on  (-L,L) • 

(iv)  For  L  >  N  a  positive  integer,  and  NL1  <  L  <  (N+1)L1#  E*(L)  consists  of  one 
isolated  positive  solution  q  and  N  j-parameter  families  Pj(L),  j  »  1,...,N,  of 
nonnegative  solutions.  If  L  *  (N+1)L1f  E  (L)  contains  one  additional  element. 

Recently  Smoller  and  tfasserman  [18]  studied  E  (L)  for  Problem  I*  in  the  case 
m  -  1.  In  contrast  to  the  result  stated  above,  they  find  that  L1  -  +•  for  m  *  1*  This 
situation  is  summarized  in  the  two  diagrams  in  Figure  1.  These  diagrams  indicate  the 
general  behaviour  of  possible  values  of  L  plotted  against  uaax  (the  maximum  of  u)  for 
u  c  E  (L) . 


Figure  1 
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A 


Having  described  the  equilibrium  set  E(L)  for  Problem  I*  we  next  turn  our  attention 
to  the  question  of  stability  of  the  various  equilibria.  As  in  the  case  m  *»  1 ,  it  turns 
out  that  both  the  trivial  solution  and  the  large  positive  solution  q  are  stable.  To 
establish  this  fact  we  apply  a  stabilization  theorem  proved  in  Part  II  in  the  general 
setting  of  Problem  I. 

We  begin  Part  II  by  proving  various  basic  existence,  uniqueness,  comparison  and 
regularity  theorems  for  Problem  I.  None  of  these  results  are  entirely  new,  but  we  know  of 
no  place  in  the  literature  where  they  are  conveniently  collected.  Further  results  of  Part 

II  -  in  particular,  the  stabilization  theorem  -  provide  us  with  a  complete  metric  space  X 

* 

of  functions  on  (0,1)  in  which  the  orbits  of  Problem  I  are  precompact.  Moreover,  if 
0  <  u0  <  1  and  u(t,uQ)  is  the  solution  of  Problem  I*  at  time  t,  then  the  urlimit 
set 

w(urt)  *  {w  e  X  :  there  exists  a  sequence  {t  },t  ♦  » 

0  n  n 

(2) 

such  that  u(t  ,u_)  -►  w  in  X} 
n  0 

is  contained  in  E(L).  For  n  »  1  this  was  proved  by  Chafee  and  Infante  [7]. 

If  E(L)  consists  of  isolated  points  only,  as  in  the  case  of  Problem  I*  and 
0  <  L  <  ,  we  obtain  that  u(t,uQ)  converges  to  a  limit  in  E(L)  as  t  ♦  •  (since 

u)(uQ)  is  connected).  If  E(b)  contains  continue  of  solutions,  as  in  the  case  for 
Problem  I*  and  L  >  ,  then  no  such  statement  has  been  proved.  However,  if  one  can  find 

a  closed  invariant  subset  K  c  x  such  that  K  n  E(L)  is  discrete,  then  for  each  u^  e  K, 

u(t,u0)  converges  to  some  point  of  K  n  £(!,)  as  t  ♦  ».  The  stability  of  the  trivial 

* 

solution  and  the  large  positive  solution  q  of  Problem  I  for  L  >  LQ  are  proved  by 

exhibiting  suitable  invariant  sets  K. 

* 

Part  I  Problem  I 

In  this  part  we  shall  consider  Problem  I,  assuming  throughout  that 
(1.1)  f (u)  »  u(1-u)(u-a),  0  <  a  <  (m+1)/(m+3) 

First  f(L)  is  studied  and  then,  calling  upon  results  from  Part  II  below,  io(uQ) 
determined  for  various  choices  of  L  and  Uq. 


Section  1  Equilibrium  Solutions 

We  are  interested  only  in  nonnegative  solutions.  A  function  v  :  [-L, L]  R+  *  [0,«) 
is  an  equilibrium  solution  of  Problem  I  when  it  is  a  solution  of  the  problem 


(ID 


{(v"V  +  f ( v)  -  0  in  <-L,L)  , 

v(±L)  -  0  . 


While  we  will  write  (II)  in  the  above  form,  v  >  0  is  called  a  solution  of  (II)  exactly 
when  w  *  v®  is  a  classical  solution  of  w*  +  f(w^m)  *  0,  w(±L)  -  0.  Clearly  v  =  0  is 
always  a  solution  of  Problem  II.  As  we  shall  show  below,  there  are  also  nontrivial 
solutions  provided  that  L  is  sufficiently  large. 

Suppose  v  ia  a  positive  solution  of  Problem  II,  i.e.  it  is  a  solution  and  v  >  0 
on  ( -L, L) .  Then  there  exists  a  £  e  (-L,L)  such  that  0  <  v(x)  <  v(^)  for  x  c  <-L,L) 
and,  clearly,  v'(£)  “  Conversely,  let  us  seek  conditions  on  £  e  (-L,L)  and  p  c  R+ 
which  guarantee  that  the  solution  of  the  initial  value  problem 


(II*) 


((v"V  +  f(v)  -  0 
v(  c)  *  u,  v’  (C)  *  o 


is  also  a  positive  solution  of  Problem  II. 

If  p  =  1  then  v  =  1  is  the  unique  solution  of  (II*)  (recall  (1.1)  -  f(1)  *  0).  If 
p  >  1,  then  f(v)  <0  for  v  >  1  implies  that  any  solution  of  (II*)  is  convex  on  its 
domain  of  definition  and  hence  cannot  satisfy  v(±L)  *  0.  Thus  (II*)  has  no  solutions 
satisfying  the  boundary  conditions  unless  p  <  1.  Consequently,  we  consider  only 
p  e  (0,1). 

To  solve  Problem  II*  we  integrate  the  equation  in  the  usual  way.  Multiply  the 
equation  by  (v*)',  and  integrate  the  result  using  the  initial  conditions  to  find 


(1.2) 


^  (vm)  '  '*  ♦  mF(  v) 


mF(  p) 


where 
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More  explicitly. 


F(v)  -  sm  ^(sjds 


F(v)  “  -  ^73  v“f16(v) 


.  2  .  .  nH-3  a  m+3 

G(v)  -  v  -  (1+a)  —  v  +  a  — 


Since  f  >  0  on  (a, 1 ),  P  is  strictly  increasing  on  (a,1).  Thus,  if  y  >  a,  we  can 


Integrate  (1.2)  to  obtain 


/?  /p 

2  J  v 


2  Jv  /F(y)-F(n) 


«Jn  -  I  C-xl 


The  integrand  in  (1.4)  has  a  singular  point  at  n  *  y,  but  F(y)  -  F(  r\)  >  6(y-r\)  for  some 

6  >  0  and  n  near  y  so  the  singularity  is  integrable.  Equation  (1.4)  defines  v 

implicitly  as  a  function  of  |i;-x|  so  long  as  v  <  y. 

If  F( y)  <  0,  then  there  exists  a  unique  v  e  (0,a)  such  that  F( v)  *  F( y)  and 

F( n)  <  F(y)  for  n  e  (v,y).  In  this  case,  (1.4)  represents  a  periodic  solution  of 
(v®)"  +  f(v)  *  0  whose  values  lie  in  ( v, y) .  Thus,  in  order  that  (1.4)  represent  a 
positive  solution  of  Problem  II  it  is  necessary  that  F(  y)  >  0. 

The  sign  of  F  is  determined  by  the  sign  of  G  and  one  checks  that  G  has  a  unique 
root  a  e  (a,1)  if  and  only  if  G(1)  <  0  or 


In  particular,  if  (1.5)  holds  then 
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F  <  0  on  ( 0 , a)  and  F  >  0  on  { a, 1 )  , 
whence  we  may  restrict  our  attention  to  the  range  a  <  y  <  1 . 

For  u  G  [a,1)  we  have  F( n)  <  F( y)  for  all  n  e  (0 ,  y).  Thus  we  can  extend  the 
integration  in  (1.4)  down  to  v  -  0.  Define 

(1*6)  Xtu)  -  /|  /o  7?(gl-F(nl  dn'  “  *  w  <  1  • 

If  y  *  a  the  integrand  in  (1.6)  may  have  a  second  singularity  at  n  *  0.  However, 

-F(n)  >  for  some  6  >  0  and  n  >  0  near  0,  so  n™"" 1  ( -F(  n)  )^2  < 

near  n  *  0.  Since  m  >  1  this  singularity  is  integrable  and  X  is  well-defined  on 
[  a,  1 ) . 

For  a  positive  solution  v  of  Problem  II,  v  *  0  only  at  ±L.  Therefore 

X(y)  -  lc-L|  *  U+L| 

from  which  we  conclude  that  £  -  0.  To  summarize,  we  have  proven  the  following  result: 
Proposition  1  Suppose  0  <  a  <  (nH-1  )/(nH-3 ) .  Then  v  is  a  positive  solution  of  Problem  II 
if  and  only  if 

/v(x)  TffarrTG  dn  -  1x1  for  1x1  <L  ' 

where  y  e  (a, 1)  and  L  e  R+  are  related  by  the  equation 
(1.7)  MU)  -  L 

and  a  is  the  unique  root  of  F  in  (a,1). 

In  view  of  (1.7),  there  is  a  positive  solution  of  Problem  II  for  a  given  interval 
( -L, L)  if  and  only  if  L  is  in  the  range  of  X,  i.e.  L  e  X([a,1)).  When  L  -  X( y)  we 
write  v(x,y)  for  the  corresponding  positive  solution.  The  multiplicity  of  these  positive 
solutions  is  the  same  as  the  multiplicity  of  the  roots  of  X(y)  -  L,  which  is  determined 
by  the  shape  of  the  graph  of  X.  Our  next  result  shows  that  the  graph  of  X( y)  always  has 
the  general  features  indicated  in  Figure  1. 
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L1  (m) 


/ m  ra(m)  n"  1 
J  2  >0  T- F(~mf  dn 


If  a  t  (0,1/2)  then  a(m)  is  defined  and  continuous  for  m  >  1.  Moreover,  as  m  -►  1 , 

F(n?m)  >  and  the  singularity  of  (n/F(a(1);1)  -  F(n;l)  1  =  (n/F(-n/1>)  1  at 

H  =  0  is  not  integrable.  It  follows  from  Fatou's  lemma  that  lim  L^Cm)  *  ».  Thus  the 

nUI 

nonexistence  of  the  small  oositive  solution  on  sufficiently  large  intervals  is  due  to  the 
nonlinearity  of  the  diffusion. 

Proof  of  Proposition  2 
Write 


MU)  -  /-  X(u)  = 


u _ n 


m- 1 


✓F(y)-F(n) 


dn 


and  use  the  change  of  variables  t  *  n/y  to  obtain  the  expression 


/J  7TO7 


m- 1 


wT 


di 


Formally  differentiating  the  integral  yields 

(1.8)  A'(u>  =  =  At u)  -  f  f?  A1  aT  . 

"  2  °  (F(y,-F(xy))3/2 

For  u  c  (a,1)  we  have  F*(u)  =  u™  V(u)  /  0  and  it  is  not  difficult  to  verify  that  the 
integral  in  (1.8)  is  convergent,  the  equation  is  valid  and,  indeed,  A*  c  C(a, 1).  If  we 
set 

9(n>  »  2mF ( n)  -  nmf(n) 
then  A'  can  also  be  written  in  the  form 


(1.9) 


A* ( u)  = 


2  u 


fU 

}0 


Q(  p)-o(  n) 


(F<  u)-F(  n) ) 


3/2 


m- 1 

n 


dn  . 
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To  study  the  behaviour  of  A(p)  as 


U  +  «/ 


A(^  =  (^o+  0  wr3rarnfdn  s 


I/u) 


I2(u) 


For  n  <  a  <  u t 


m-  1  m-  1 

n _  *  _n _ 

7ft m)-f J n)  T^fThT 

and  we  have  already  noted  that  the  right-hand  side  is  integrable  on  [0,a],  so 

n- 1 

Vu)  *  /o  T^fTnf  dn  =  A(a)  as  11  +  0  • 

On  the  other  hand,  by  previous  remarks,  on  each  compact  subset  of  a  <  n  <;  \i  <  1 , 
nm  V/F(  u)“F(  n)  is  dominated  by  a  multiple  of  (p-n)  ^  ,  so  1  (\x)  +  0  as  y  *■ 
it  follows  that  A  e  C([a,1))« 

Since  0(a)  =  -amf(a)  <0  we  can  choose  y  >  0  such  that  9(jj)  <  0(a)/2  <  o 

U  c  [ a,  a+-y).  In  addition,  since  9(0)  =0  we  can  choose  5  >  0  such  that 

I  9(  n)  I  <  -0(a)/4  for  n  e  (0,5).  Thus,  in  particular,  \i  e  [a,  a+y)  and  n  e  [0, 

imply 

e<g)  -  e(n)  <  <  o  . 

4 


Write 


,  1  rf 5  .  rg,  9<y)-9(n)  n*-1  , 

AMu)  =  ^  {/0  +  r5>  — — — 372  n  dn  =  +  J2(p)  • 

(  f( -f( n) ) 


Arguments  like  those  above  show  that  remains  bounded  as  u  4-  a.  For  u  f  (a 


we  have 


■V-  <  ^  /.* 


<F( y )  -F(  n) ) 


*  in- 1  -  3/2 

As  u  +  a,  the  integrand  of  (J  (y)  converges  pointwise  to  rj  {-F(n))  wh:ch 

.  .  ...  m- 1  -  3/2  (m+1)  -1/2  (m+r;)  .  .  .  .  , 

behaves  like  n  ( n)  =  r\  near  n  =  0  and  so  is  not  mtegrable.  Thus 

* 

Fatou's  lemma  yields  J  (y)  -*-+<»  (and  so  A' (  y)  ♦  -°°)  as  p  *  a. 

To  show  that  A(y)  ♦  »  and  A1  (  y)  *  00  as  u  t  1  we  argue  similarly.  First,  the 
integrand  in  A(  y)  tends  to  n™  1  //?  ( 1  ) -F  (  n)  •  Now  F  *  (  1  )  *  f  { 1  )  =  0  and  F"  {  1  )  <  0,  so 

this  has  a  nonintegrable  singularity  at  1  and  it  follows  that  A(y)  ♦  +  «’  as  y  t  1.  Th< 
inteqrand  in  A*(y)  also  tends  to  a  nonintegrable  limit  since  9'(1)  >  0,  and  it  follows 
that  A'  (  y)  +  +°°  as  y  +  1. 

We  now  turn  to  the  proof  of  (iv).  For  this  we  follow  closely  the  proof  given  in  the 
case  m  *  1  by  Stttoller  and  Wasserman  [18]. 

To  begin  with,  we  note  several  properties  of  the  function  6.  They  are  all  proven  bv 
elementary  calculations  using  the  explicit  formulas  for  f  and  F.  We  shall  omit  the 
details . 

A.  There  is  a  y^  e  (a,1)  such  that  0( n)  <  0  on  (0,y^) 
and  Q(  n)  >0  on  (  y^  ,  1  ]  . 

B.  There  is  a  ^  e  (0,^)  such  that  O'(n)  <  9  on  (0,y^) 
and  0* ( n)  >  0  on  ( , 1 ] • 

C.  There  exists  a  y^  e  (0,y^)  such  that  (  n0*  (  n) )  *  <0  on  ( 0 , y^ ) 
an  f  n) ) *  >  0  on  ( p^ , 1 1 . 

It  follows  from  properties  A  and  B  that 

A'(y)  >  n  on  (u  ,1) 

and,  if  a  <  y^' 

A'  (  y )  <  0  on  (  n,  y^  ) 

Thus  we  need  only  consider  A*  for  max(a,  y^ )  <  u  *•  • 

To  proceed  we  need  to  examine  AM*  For  this  purpose,  let 
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o  <,  n  <  u 


and 


(  6jh)  (  n)  *  h(  U)  •  h(  n)  , 

(<5  h)(n)  *  yh(  p)  -  nh(  n) »  o  <  n  <  p 


(«3h)(n)  =  u"h<u>  -  nmh(n),  o  <  n<  u  . 


A  computation  similar  to  the  one  yielding  (1.9)  produces 


.  (6  F>(6  0'  >-  -  <6  0)<6  f) 

..  I  fU  I  2  21  3  m- i  _  m- 1  ...  . 

(1.10)  A"(  p)  =  — j  Jo  - X/7 -  ^  +  k  (  V) 

2\i 


(^F) 


Adding  (K  -  — jp)A*(u)  to  both  sides  of  (1.10)  we  obtain 


,  .  (  6  F)  (  6  0* )  +  6  etKp&F  -  -  S,f>  . 

.......  m-1JJt,  .  1  rU  1  2  1  1  23  m- 1  . 

C  u>  +  (K - ”)  A'(p)  =  — -  /0  - 575 -  0  dr) 


2p 


(  6^) 


Observe  next  that  6^0  =  2m6^F  -  6^f  so  that  KpS^F  -  |  S^f  =  (Kp  -  301)6^  +  |*  5^.  Thus, 
if  we  choose  K  -  3m/ p  we  arrive  at  the  expression 


,  |  («,  0)2  +  (S.FKS.e1  )  m  , 

v  .  2m+  f  . , ,  v  1  rV  2  1  1  2  m-1  „ 

- - - a. 


In  view  of  properties  B  and  C  above,  we  have 

620*  ■  u0'(g)  -  n9'  (  n>  >  0  for  max(.-,p2>  <  p  <;  1 
and  it  follows  that 


(1.11) 


AM<  p) 


2m+1 

U 


AMu)  >  0 


for  max(a,  u2>  <;  u  < 


As  noted  above,  A?(u)  <  0  for  p  <  max(a,  p^)  and 
Therefore  A'  has  at  least  one  zero  in  the  interval  J  = 
(1.11)  implies  A*(p)  >  0  at  any  such  zero  and  therefore 
completing  the  proof. 


A1  (  p)  >  0  for  p  >  . 

lmax(a,p2),  u  ].  The  relation 
there  can  be  at  most  one, 
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Remark:  As  the  proof  shows,  the  result  (iv)  is  dependent  only  on  the  properties  A,  B,  C 
of  0. 

Propositions  1  and  2  provide  a  complete  characterisation  of  the  set  of  positive 
solutions  of  Problem  II.  For  LQ  <  L  let  y+(L)  denote  the  largest  solution  of 
L  -  X(y)  and  for  LQ  <  L  <  let  y  (L)  be  the  smallest  solution  (so 
y+(LQ)  *  MJV>.  We  <3ist^n^u^sh  the  following  cases 

0  <  L  <  Lq.  There  are  no  positive  solutions. 

L  =  Lq.  There  is  a  unique  positive  solution  v(  •  ,  y+(L^  ) )  • 

LQ  <  L  <  L^.  There  are  two  positive  solutions  p(  *,L)  -  v(*,y  (L)) 

and  q( • ,  L)  »  v(*,y+(L))  with  p  <  q  everywhere  on  (-L,L). 

L  >  Lj.  There  is  one  positive  solution  q(  *,L)  *  v(«,y+(L)). 

Since  v(*,y)  depends  continuously  on  y  and  y±(L)  are  continuous  on  their  domains, 
p  and  q  are  continuous  functions  of  L  on  their  domains. 


We  now  show  that  v(*,a)  -  v(»#U_(I^))  generates  families  of  nonnegative  solutions  of 
Problem  II  on  intervals  ( -L, L)  with  L  >  .  For  y  c  (a,  1]  we  have  F(y)  >  0  so  that, 

according  to  (1.2)  ( vm) ' ( ±X( y ) , y )  ?  0.  However,  F(a)  “  0,  so  ( vm) 1 (±X( a) , a)  * 

(vm)'(±L1,a)  *  0.  It  follows  that  v(x,a)  extended  as  0  for  L  >  |x|  >  is  a 

solution  of  Problem  II  for  L  >  L^  and  so  is 


r(x;h)  « 


{v(x-h;a)  for  |x-h|  <  L^ 
0  for  | x-h  j  >  L1 


provided  |h|  <  L  -  L  •  More  generally,  we  may  piece  several  such  solutions  together  if 
their  supports  are  disjoint.  Let  N  be  a  positive  integer  and  L  >  NL^.  For  each  N- 
vector  c  -  which  satisfies 


the  function 


r(x»£) 


v(x-£i#a)  for  |x-q|  <  L<|  , 


if  lx-^1  >  L1  for  i  -  1,...,N 


la  a  nonnegative  solution  of  Problem  II.  We  shall  use  PN(L)  to  denote  the  collection  of 
functions  v(«,£)  where  £  e  RN  satisfies  (1.12). 

Clearly  a  nonnegative  solution  of  Problem  II  is  either  positive  or  belongs  to  some 
Fn(L).  We  thus  have: 

Proposition  3  For  L  >  let  N  be  the  integral  part  [ L/L ^  ]  of  L/L^.  Then  with 
N 

P (L)  *  U  P  (L)  we  have 
j-i  3 

E  (L)  -  {q(  .,L)  }  U  P(L)  • 

Remark.  If  L/L 1  >  [L/L^]  •  N  then  PN(L)  is  a  true  N-parameter  family,  while  if 
N  -  L/Ly ,  PN(b)  contains  only  r(x?$)#  ^  =  (i-DL^  +  1^/2.  Combining  Propositions  1.1, 
1.2  and  1.3  we  arrive  at  the  complete  description  of  E(L)  given  below: 

Theorem  4. 


for  0  <  L  <  Lq  , 


E(L)  *  <  {0,  p(  •;!.),  q(»|L)}  for  LQ  <  L  <  , 

I  (0,  q(«,L)}  <j  P(L)  for  <  L 


Section  2.  Stability  Theory 

We  now  turn  to  the  question  of  the  large  time  behaviour  of  the  solution  of  the 
initial-bounday  value  problem 


r  u.  *  (u  )  + 

[  t  XX 

f(u)  in 

( -L,L)  x 

<I#) 

\  u(±L,t)  *  0 

in 

+ 

R 

\  u( *,0 )  -  uQ 

in 

(-L,L) 

where  f 

is  given  by  (1.1)  and 

(1.13) 

00 

Uq  e  L  (“Ii,L)  f 

0  «  u0  < 

1  a .  e , 
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In  what  follows  we  shall,  for  convenience,  write 


$1  -  (-1*,  L)  ,  *  n  X  (0,T)  ,  Q  -  «  x  R+  . 

Definition*  A  solution  u  of  Problem  I*  on  [0,«0  is  a  function  u  :  [0,*)  ♦  L*(fl) 
with  the  properties 

(i)  u  eC([0,»)  :  L^ft))  n  l"(Qt)  for  T  >  0. 

(ii)  J  u(T)*(T>  -  //  (u*  +  uV  )  -  L  u  *(0)  +  //  f(u)*> 

H  Z  XX  U  U 

2  — 

for  all  T  >  0  and  C  (Q)  such  that  yP  >  0  in  Q  and  \p  -  0  at  x  -  ±L. 

A  subsolution  (supersolution)  of  Problem  I*  on  (0,®)  is  a  function  satisfying  (i) 
and  (ii)  with  equality  replaced  by  <  (respectively,  >). 

Theorem  5  (Existence  and  comparison) 

* 

(i)  If  (1.13)  holds  Problem  I  has  a  unique  solution  u  on  [0,») 
and  0  <  u  *  1  a.e. 

(ii)  If  u  is  a  subsolution  and  u  is  a  supersolution  of  Problem  I* 
then  u  <  u  a.e*  on  Q. 

Theorem  5  is  a  consequence  of  more  general  results  in  Part  II  (Theorem  12  and  13).  We 
will  denote  the  solution  of  Problem  I*  with  the  initial  value  uQ  by  u(t,u0).  Let  X  be 
the  complete  metric  space  with  the  metric  d  given  by 


(X  -  {u  e  L*(n)  :  0  <  u  <  1  a.e.  and  um  e  H^(fl)}  , 

. 

d(u,v)  »  lu-vl  +  Mu™1  «  vm)  I 

L  (ft)  X  L2 ( fl) 

Recall  the  definition  of  u>(u0)  (Introduction,  eg.  (2)).  A  stablization  theorem  (Theorem 
18)  is  proved  in  Part  II  which  applies  to  Problem  I*  to  yield: 

Theorem  6  (Stabilization):  Let  Uq  satisfy  (1.13).  Then  (u(t,u^)  :  t  >  1}  is  a  compact 
subset  of  X,  w(uQ)  is  nonempty  and  connected  in  X  and  oj(Uq)  c  ^(L). 

We  will  use  these  results  to  show  the  stability  of  the  equilibrium  solutions  v  *»  0 
and  v  -  q.  For  this  purpose  we  introduce  the  notion  of  sub-and  supersolutions  of  Problem 
II# 

(vm)"  +  f ( v )  -  0  in  ( -L,L) ,  v(±L)  *  0  . 
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A  (weak)  subsolution  of  Problem  II  is  a  function  v  e  C([-L,L])  for  which 

+  ^f(v))dx  >  0  for  all  *  e  C2(ft),  ^  >  0  and  *  (±L)  »  0  and  v(iL)  <0.  A 
(weak)  supersolution  is  defined  by  reversing  the  inequality  and  requiring  v(±L)  >  0. 

Let  v  and  v  be  respectively  a  sub-and  a  supersolution  of  Problem  II  and  let 

—  00  — 

[v,v]  =  {w  t  L  (ft)  :  v  <  w  *  v  a.e.  on  ft) 

Proposition  7  Let  uQ  e  [v,v]  satisfy  (1.13).  Then 

(i)  u(t,uQ)  e  [v,v]  for  t  >  0 

and 

(ii)  ortUg)  c  [v,v]  n  x* 

Proof.  It  follows  from  the  definitions  that  v(v)  is  a  time-independent  subsoj.ction 
(supersolution)  of  Problem  I*.  Hence  Theorems  5  and  6  imply  (i).  The  assertion  (ii) 
follows  from  Theorem  6  and  the  fact  that  ly,v]  n  X  is  closed  in  X. 

Corollary  8.  If  uQ  e  (v,v)  satisfies  (1.13)  and  [v,v]  n  f  =  (g)  is  a  singleton,  then 
u(t,Ug)  >  g  in  X  as  t  ♦ 


We  next  give  three  applications  of  this  Corollary  to  the  determination  of  domains  of 
attraction  of  the  various  isolated  elements  of  E(L). 

(1.15)  Let  L  e  (L  ,L  1.  Choose  £  e  [L  ,L)  and  5  e  (-L,L)  such  that  ~L  ^  -  i, 

£  +  £  <  L.  Set 


v  ( x )  = 


{p  (  x-  £ ,  £ )  for  X  e  l  £-£,£+£)  , 

0  for  x  l  fc-£,£+£l 


It  now  follows  from  Corollary  8  that  u(t,uQ)  +  q(L)  in  X  as  t  ♦  «  whenever 
u0  €  f£,vj. 
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(1.16)  Let  Le  [L0 , L1 ) .  Chose  A  e  (L,L  and  n  t  ( -L,L)  such  that 
n  -  A  <  -L,  L  <  n  +  A,  and  let 

w(  x)  -  p( x- n,  A) 

Then  w  is  a  supersolution  of  Problem  II.  Since  »  0  is  a  subsolution# 
invariant.  In  this  case  (see  Figure  3(b)) 

w(  0 )  *  p(-n,A)  <  P<0,A)  <  p(0,L>  . 


Hence 

[w,w]  n  E(l)  -  {0} 

and  therefore  u(t,Ug)  +  0  as  t  +  •  in  X  when  u^  e  [w,w]. 
(1.17)  Let  L  >  Lr  For  x1  ,x2  e  ( -L , L )  such  that 

-L  <  x,  -  L,  <  x2  -  Ll  <  x,  +  t.,  <  x2  +  <  L 

define 


v+(x)  *  max(p(x-x1  ,L1  ),  pfx-x^L^) 
v  (x)  **  min(p(x-x1  ,L<J  ) ,  ptx-x^L^) 


See  Figure  4.  It  is  easy  to  verify  that  v+  is  a  subsolution  and  v~  is 
of  Problem  II.  Moreover, 

[v+,1]  n  E(L)  =  (q(L)},  [0,v”)  n  E (L)  -  {0}  . 

Thus 


lim  u(t,uQ)  *  q(L)  for  all  uQ  e  [v*,1] 

t-H» 


and 


lim  u(t,uQ)  -  0  for  all  uQ  e  (0,v  ] 
t-H« 


[w,  w]  is 


supersolution 
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It  is  clear  that  more  elaborate  domains  of  attraction  for  q  and  0  can  be 


constructed.  We  leave  further  constructions  to  the  interested  reader. 

Part  II.  General  Theory 

In  this  part  of  the  paper  we  first  prove  existence,  uniqueness  and  continuous 
dependence  on  initial  data  of  solutions  of  Problem  I  of  the  introduction  together  with 
comparison  results.  These  results  are  all  more  or  less  known  in  various  contexts,  but  the 
presentation  here  collects  them  quite  conveniently.  (See  the  remarks  at  the  end  of  this 
section.)  After  this  the  stabilization  result  used  in  Part  I  is  proved. 


Section  3.  A  Preliminary  Estimate 
Consider  the  problem 


(III) 


u  =“  n(u)  +  g(x,t) 
t  xx 

\  u(±L,t)  *  ^±(t) 

V  u( x,0 )  -  uQ(x) 


(x,t)  e  Q  , 
t  e  (0, »)  , 

x  e 


where  we  assume  that  n  and  the  data  g,  \J/±  and  Uq  satisfy  the  following  set  of 
assumptions : 

A1 .  n  s  R  ♦  R  is  locally  Lipschitz  continuous  and  nondecreasing, 

A2.  g  e  l\qt)  for  each  T  >  0, 

A3.  4>±  e  LJoc(  >' 

A4.  uQ  e  L“(n). 

These  will  be  called  "assumption  A". 

Definition.  A  solution  u  of  Problem  III  on  [0,T)  is  a  function  u  with  the  following 
properties : 


(i)  u  E  C(  [0,T1  :  i/fnn 


n  L  (Ql 
T 
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(ii)  J  u(t)tft)  -  jj  (utf  +  n(u)^  )  + 

“  Q  t  XX 

*t 

(l  n( Cs )  )^x(L,s)  -  n(<Ms)  s)ds 

3  /o  v(0)  ♦  // 

Qt 

2  — 

for  all  ^  e  C  (Q  )  such  that  ^>0,^=0  at  x  *  ±L  and  0  <  t  <  T.  A  solution  on 
T 

(0,<»)  means  a  solution  on  each  [0,T] ,  and  a  subsolution  (supersolution)  is  defined  by 
(i),  and  (ii)  with  equality  replaced  by  <  (>). 

Proposition  9  Let  u  be  a  supersolution  of  Problem  III  on  [0,T]  with  data  g,  uQ,  ^ 
and  u  be  a  subsolution  on  [0,T]  with  data  g,  uQ,  ^  all  satisfying  assumption  A. 
Then  if  we  have  for  each  X  >  0  and  0  <  t  <  T 


(2.1  ) 


Xt 


(u(t)  -  <  /  (u0-u0>+  ♦  //  e*S(g-g  +  X(u-u))  + 


where  r+  -  max(r,0). 

Proof.  Since  u  (u)  is  a  supersolution  (subsolution)  and  if  >  0  (<  0)  at  x=  -L(+L), 

A 

we  find,  using  n(  4^. )  >  h(i|^)  (by  Al) 


(2.2) 


(u(t)  -  u(t)K(t)  -  JJ  (u  -  u)(^fc  ♦  aV?xxJ  < 

«  fa  (u0  "  +  //  ( g  -  §V  > 


where  a  =  (n(u)  -  r)(u))/(u  -  u)  for  u  +  u  and  a  =  0  otherwise,  for  all  ^  e  C2(Qt) 
such  that  ^  >  0  and  =  0  at  x  *  ±L  and  0  <,  t  <  T.  By  Al  and  the  boundedness  of 

/>  CD 

u,  u,  we  have  a  e  L  (Q^)  and  a  >  0. 

We  now  construct  a  special  sequence  of  functions  }  to  use  in  (2.1).  Fix  T  >  0 

n 

and  choose  a  sequence  {a  }  of  smooth  functions  such  that 

n 


and 


<  a 


n  n 


<  Kal 

00 

L  <qt> 


+ 


n 
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2 

(a  -a)//a~  >0  in  L  (O  ) 
n  n  T 

This  is  easily  seen  to  be  possible.  Next  let  x  e  C~(ft)  ***  such  that  0  <  X  <  1*  Finally 
let  be  the  solution  of  the  backward  problem 


f  ^  .  +  a  v?  * 

nt  n  nxx  n 

for 

X 

e 

fl,  t  e  (0,t) 

(2.3)  ^ 

*  (±L,t)  -  0 
n 

for 

t 

e 

10, T)  , 

l 

k  *>n(x,T)  -  x<x> 

for 

X 

e 

n  . 

This  is  a  nondegenerate  parabolic  problem  and  has  a  unique  solution  V  n  e  c“(QT>. 

Lemma  10  The  function  has  the  following  properties: 

(i)  0  <  <  •X(t-T)  on  Q  , 

n  T 

(il)  //  an^nxx>2  <  C  ' 

CT 

(iii)  sup  L  )2(t)  <  c 

a  »*  nx 

0  <t  <T 

where  the  constant  c  depends  only  on  x* 

Proof.  Part  (i)  is  immediate  from  the  maximum  principle  and  0  <  x  <  1*  To  prove  (ii)  and 
(iii),  multiply  the  equation  solved  by  by  ^nxx  and  integrate  over  fi  x  (t,T)  to 

find  -  after  an  integration  by  parts  - 


-£  4  Vnxt  +  £  /ft  VVJ  ’  *X  /ft 


I  /ft  (%x,2,t)  +  /t  /ft  VO"  +  X  /I  /ft  (%x>2  ■  I  /ft  %> 


from  which  we  have  the  desired  estimates. 

If  we  set  t  -  T  and  in  (2.2)  we  obtain 


1 


£ 
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(2.4) 


L  (u(T)  -  u(T))x  -  //  (u-u)  (a -a  )^  < 

J  il  n  n  nxx 

<  (uQ  -  ^0)^n(°)  +  //  (g-g  +  A(u-u))^n 

PT 

<  J  (u0  -  a0)+  e'XT  +  //  ( g-g  +  X(u-u))  +  eMs‘T) 
«T 

Since 


// 


|a-a  I 
n 


// 

qt 


|a-an! 


( /a”  1^? 
n 


nxx 


I) 


we  have ,  by  Lemma  10  ( i i ) , 


(a-a  ) 

l(a-a  )*>  I  .  <  I—7  n  I  ,  t/a"  *  I  , 
n  nxx  1  /a  ,2  n  nxx  2 

L  n  L  L 


v,  <•-•-> 


/a  2 
n  L 


which  tends  to  zero  as  n  «•  by  the  choice  of  an.  (The  spaces  L1  and  L2  here  are 
taken  over  ftp.)  Thus,  letting  n  ♦  ®  in  (2.4)  we  obtain 


(2.5)  fa  (u(T)  -  u(T) )  x  <  f0  (uQ  -  uQ)+  e_XT  +  //  (g-g  +  X(u-u))  +  cX(s*T)  . 

qt 

This  inequality  holds  for  every  x  e  CQ(ft)  with  0  <  x  *  1*  Hence  it  continues  to  hold 
for  x<x)  *  1  on  (x  •  u(T)  >  u(T) }  and  x  =  0  otherwise  (i.e.  x  =  sign  (u(T)-u(T) ) +  , 
completing  the  proof.  (Clearly  T  may  be  replaced  by  any  t,  0  <  t  C  T  in  the  argument.) 
Corollary  11.  Let  u  and  u  be  solutions  of  Problem  III  with  data  g,  uQ,  ^  and 
g,  Uq  ,  Then 

(2.6)  ,luft)  -  u(t)l  <  !u  -  u  I  +  J*  lg(  s )  -  g(s)l  ds  . 

L  { fl)  L  (ft)  L  (ft) 
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Thus,  in  particular,  solutions  of  Problem  III  are  unique. 
Proof.  Setting  X  =  0  in  (2.1)  we  obtain 


fa  '“(t)  -  “<t))+  <  fa  (uQ  -  uQ)  +  +  //  <g  -  g 


V 


and 


/fl  t fc>  -  u(t))+  <  /n  (u0  -  u0)+  +  //  (g  -  g)'1 


Adding  these  estimates  yields  (2.6). 

Uniqueness  and  continuous  dependence  for  Problem  II . 


m-1 


We  now  return  to  Problem  I.  Here  we  take  n(u)  ■  |u|  u  which  satisfies 
solution  of  Problem  I  we  mean  a  solution  u  of  Problem  III  with  g  «  f(u),  ■ 

on  for  sub-and  supersolutions.  Since  solutions  are  bounded  by  definition , 

oo  1 

g  e  L  (Q^)  c  l  (Q^),  and  we  may  use  the  previous  results. 

Theorem  1 2 

(i)  Let  u,  u  be  solutions  of  Problem  I  on  [0,T]  with  initial  data  uQ 

and  Uq  respectively.  Let  K  be  a  Lipschitz  constant  for  f  on  [ -M,  M] 

where  M  =*  max(lul  r  lul  ).  Then 

L*{0t) 

I  u  ( t )  -  u(t)l  <  eKt  Bu  -  u  B 

tVfl)  °  °l\«) 


(ii)  Let  u  be  a  subsolution  and  u  a  super solution  of  Problem  I  with  initial 
Uq  and  Uq.  Then  if  u^  <  u^  i:  follows  that 

.  A 

u  <  u 

Proof.  With  the  assumptions  of  (ii).  Proposition  7  yields 


(2.7)  e 


Xt  r 


'  ^  (u(t)  -  u(t))  <  (uQ  -  uQ)  +  /q  e^S(f(u)  -  f  ( u)  +  X(u  -  u] 


A1 •  By  a 
0  and  so 


data 


+ 

) 
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Set  X  *  K  (defined  as  m  (i)).  Then  r  >  f(r)  +  Kr  is  nondecreasing  on  [ -M , M]  and 


Thus  if  we  write 


(f(u)  -  f(u)  +  K(u  -  u))  <  2K(u  -  u) 


h(t)  =  eKt  fQ  (u(t)  -  G(t))+ 


(2,7)  implies 


h(t)  <  h( 0 )  +  2K  J*  h(s)ds 

2Kt 

which  implies,  by  Gronwall's  lemma,  that  h(t)  <  h(0)e  or 

Jfl  <“<*>  -  “(t))+  <  eKt  fa  tuQ  -  uQ>+  . 

This  proves  (ii).  The  assertion  of  (i)  follows  by  adding  the  corresponding  inequality  for 
(u  -  u)+. 

Remark.  The  proofs  here  are  correct  for  any  q  satisfying  A1  and  do  not  require 

0  <  a0,  G0  <  1. 

Section  4  Existence 

We  begin  by  regularizing  the  problem.  Let  e  >  0  and  consider 

f  u  *  (um)  +  f  (u)  ( x , t )  e  Q_ 

t  xx  e  T 

u(±L,t)  -  e  t  e  (0 ,T] 

s.  u( x , 0 )  =  uQ(x)  +  e  x  e  ft 

f  (u)  -  f(u  -  e)  • 
e 

The  properties  of  f  and  uQ  we  will  use  are: 


(I  ) 
e 


where 


r 


(H)  f  :  R  +  R  is  locally  Lipschitz  continuous,  f(0)  -  f (  1 )  -  0 
and  e  0  <  uQ  <  1, 

which  we  refer  to  as  hypotheses  H.  For  a  while  we  also  assume  uQ  e  Then 

2i  — 

Problem  I  has,  by  classical  results,  a  unique  solution  u  e  C  '  (Q  )  and 
£  £  T 

(2.8)  e  <  u  <  1  +  e  in  Qm  . 

£  T 

Multiplying  the  equation  of  1^  by  {u™)^  and  performing  obvious  manipulations  yields 


ro+1  2 


4m  rt 


f\  fa  (u  *  2  fa  K>x“»‘  -  m  fa  FE(uE)(t> 


(m-M  ) 


(2.9) 


2  fa  {(ue)x(T)}"  -  m  fa  FE(uE)(T) 


for  0  <  t  <  t  <  T,  where 


F  (u)  -  f1  Sn_1f  (s)ds  . 
£  J  t  £ 


In  particular,  putting  r  m  0, 


(2.10) 


sup  /  {(u1")  (t)}2  <  K 

0  <t  <T  M  C 


where  K  depends  on  f. 


r  ,  ®.2 

Vuo\ 


but  not  on  e  e  (0,1)  or  T.  Set 


Then  (2.9), 


(2.10)  imply 


(2-11) 


e <  V  <  ( 1  +  e) 

E 

iv  (t) a  _  <  /k 

05  L2<«) 


IV 


et 


t,2(0,T:L2(n) 


d  U  l«” 


m+1  2m 


m+ 1 


K  . 


2 

It  follows  from  (2.11)  that  {v^},  0  <  E  <  1  is  equicontinuous  from  [0,T}  into  L  (ft) 

1  2 

with  values  in  a  bounded  subset  of  H  ( ft)  (which  is  compactly  imbedded  in  L  (ft)).  Hence, 

2 

by  Arzela-Ascoli' s  Theorem,  there  is  a  v  e  C([0,TJ:  L  (ft))  and  e  +0  such  that 

n 

v  v  in  C(  (0,T]  :  L2(ft)).  Then  u  •*>  =u  and  um  ♦  um  in  C(  [0,T]  :  L2(ft)). 

e  EE 

n  n  n 

It  is  very  simple  to  show  that  u  is  a  solution  of  Problem  I,  and  we  omit  this.  (Note 

that  C([0,T]  :  L2(fl))  c  C(  fO.T)  :  thfl)).) 

It  remains  to  remove  the  restriction  uQ  e  C*{ft).  To  this  end,  let  uQ  satisfy  (2.8) 

and  choose  a  sequence  {\in  }  c  C*(ft),  0  <  u.  <1,  such  that 

On  0  On 

(2.12)  luft  -  u_  I  «  ♦  0  as  n  ♦  »  . 

0  0n 

Let  be  the  solution  of  Problem  I  with  initial  data  Upn.  According  to  Theorem  12(i) 

KT 

(2.13)  sup  lu  (t)  -  u  (t)l  <  e  Hu  -  u  II 

0 <t <T  3  L  (ft)  3  L  (ft) 

where  K  is  now  a  Lipschitz  constant  for  f  on  fO,T]  .  By  0  <  <  1,  (2.12),  (2.13) 

there  is  a  u  e  C( [0,T]  :  L1 (ft) ) ,  0  <  u  <  1,  such  that  u^  +  u  in  C( (0,T]  :  L1  (ft)  ) . 
Clearly  u  is  a  solution  of  Problem  I  and  so  we  have  proved: 
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Theorem  13*  Let  hypotheses  H  hold.  ihen  Problem  I  has  a  unique  solution  u.  Moreover, 


0  <  u  <  1. 


Section  5.  Regularization 

Xn  this  section  we  shall  prove  a  regularizing  property  of  the  solution  operator  of 
Problem  I. 

Theorem  14.  Let  hypotheses  H  hold  and  u  be  the  solution  of  Problem  I.  Then  for  each 

t  >  0  there  is  a  constant  independent  of  Uq,  such  that 

(i)  um(t)  e  L  (ft)  for  t  >  t 
x 

and 


flu"(t>  I  <  M  and  essential  variation  um(t)  <  M 

T  x  T 

for  t  >  x 

Proof .  Following  [6]  we  denote  the  solution  operator  of  Problem  III  with  n(r)  *  |rj  r 
and  ^  *  0  by  S(t,uQ,g)  -  that  is,  S(t,Ug,g)  is  the  solution  of  Problem  III  at  time  t 

if  this  problem  has  a  solution.  By  Corollary  11,  S  has  the  Properties 


MS(t,uQ,g)  -  S{ t ,uQ 


50I  1  ♦  J*  «g(s)  -  g(s) i  ,d. 


and 


(2.15) 


1  1  m 

Xm  1  S(Xt,uQ,g)  *  S(t,  X™""*  uQ,  X*-1  g^) ,  X  >  0 


where  g  (t)(*)  *  g(Xt)(«).  To  establish  (2.15)  one  merely  checks  -  using  n(r)  -  |r|m  'r 
7  1  1 

-  that  Xm  1  S(Xt,uQ,g)  is  indeed  a  solution  of  problem  III  with  data  Xm-1uQ,  X™" 1 g ^  in 

place  of  uQ,g  ar*d  1‘nvokes  fc^e  uniqueness.  Now  the  solution  u(t,uQ)  of  Problem  I  is 
exactly  a  solution  of  u  *  S( t ,Uq , f (u) ) .  By  Theorem  7  of  [6],  properties  (2.14)  and  (2.15) 
of  S  and  the  Lipschitz  continuity  of  f  imply  that  for  t>0,  0<h<T,  t>0 


£  lu(t+T+h,UQ)  -  u(t+T,uQ)l  1  »  1  lu< i+h,u( t,uQ) >  -  u( T,u(t,UQ) ) I  ^ 

L  L 

<  \  E(t,  lu(t,u0)l  ,) 

L 

where  E  is  a  nondecreasing  function  of  its  arguments.  Since  0  <  u  <  1,  we  have 

lu(t,uQ)l  <  meas  ft  -  2L  and  it  follows  that  t  1E( t,2L)  is  a  Lipschitz  constant  for 
1/ 

t  *  u(t+T,uQ)  on  [0,®). 

The  proof  is  completed  by  means  of  the  follows  lemma: 

Lemma  15.  Let  v(t)  be  Lipschitz  continuous  with  constant  L,  and  v(t),  z(t)  be 
continuous  from  [0,®)  into  L^fft)  and 

v  *  w  +  z  in  V* ( ft) 

t  XX 

Then  w(t)  £  L*(Q)  for  each  t  and 
x 

(2.16)  essential  variation  w(t)  <  L  ♦  lz(t)l 

X  L 

We  apply  this  lemma  to  the  equation 

ut  -  (u”>xx  +  E<u> 

which  holds  in  V *  (i.e.,  in  the  sense  of  distributions).  As  shown  above,  t  ♦  u(t,uQ) 

is  Lipschitz  continuous  from  (t,®)  into  L*(fl)  with  a  constant  L^  independent  of  uQ. 

By  Lemma  15  um(t)  e  L°°(ft)  for  t  >  t  >  0  and  the  variation  of  um(t)„  is  bounded  by 
x  x 

L  +  lf(u( t) )  I  4,  which  is  bounded.  If  v  :  [0,1]  +  R  is  smooth  then  v  (a)  * 

T  L1 

v(1)  -  v(0)  for  some  a  e  (0,1)  and  then 

Iv  I  <  |v  (a)}  +  variation  v 

x  «  X  X 

L 

<  2lvl  +  variation  v  • 
l" 

By  approximation  with  smooth  functions  we  conclude 

l(um)  I  <  2lum|  +  ess  variation  (um) 

X  ®  ®  X 

L  L 

<  2  +  ess  variation  (um)x  • 


Thus  the  assertions  of  the  theorem  are  established 


It  remains  to  prove  the  lemma. 
Proof  of  Lemma  15.  Define  the  averages 


v 


h 


1  rt+h 
h  Jh 


v, 


i  ft+h 
h  ^t 


w, 


2  rt+h 

h  't 


z 


Indeed,  these  averages  can  be  defined  for  an  arbitrary  distribution  F  on  Q  by 

V*  - F[  i  /tt-h)^  for  * e  vQ) 

and  then  it  is  easily  checked  that  the  operation  commutes  with  differentiations.  Hence 

h'^vft)  -  v(t-h)  *  vht  -  whxx  +  2h*  We  conclude  whxx  e  I*1  and 

variation  w.  *  lw  I  ,  <  L  +  lz_  (t)| 

hx  hxx  1  h  i 

L  L 

♦  L  ♦  fz(t)  f  as  h  ♦  0  . 

L 

1  00 

Since  w^,  w^xx  remain  bounded  in  L  ( ft) ,  w^x  is  bounded  in  L  (i))  and  w^x  in 

LP(fl),  1  <  p  <  ».  Letting  h  «►  0  and  using  the  lower  semicontinuity  of  the  variation  we 
obtain  (2.16). 


Section  6.  Stabilization 


Let  0  <  uQ  <  1  and  u  «  u(t,UQ)  be  the  solution  of  Problem  I  emanating  from  u 
For  each  t  >  0  define  the  semiorbit 


0* 


Yt<V  *  (u<t,u0>  :  t  >  t}  . 

According  to  Theorem  14,  y  (u„)  c  x  where  X  is  the  complete  metric  space  consisting 

T  0  T  T 

of  those  w  e  L*(ft)  such  that 


(2.17) 


0  <  v  <  1,  (Wm)  e  lw")  I  <  M 

X  X  «  T 

L 


and  essential  variation  (w  )  <  M  , 

x  T 


-28- 


where  is  as  in  Theorem  14,  equipped  with  the  metric 

(2.18)  d(u,v)  *  lu-vl  +  «<u®  -  vm)  I 

L  («)  X  L  (0) 

One  easily  checks  that  is  complete.  Moreover,  is  compact,  Indeed, 

ffl  1  00  1 

(w  :  w  e  X^}  is  bounded  in,  e.g.,  W  ( fl)  and  is  thus  precompact  in  L  (fl).  It 
follows  that  {w  :  w  e  X^}  is  precompact  in  Similarly, 

{(wm)x  :  w  6  X^}  is  bounded  in  L~(ft)  and  in  variation.  Thus  it  is  precompact  in 

1  •  p 

L  (0)  and  then,  by  the  L  boundedness,  in  every  L  (0) #  1  <  p  <  «».  The  compactness  of 

X  follows.  We  also  let 
T 

(2.19)  X  -  {u  e  L°°(n)  :  0  <  u  <  1,  <utt)x  e  L2(«)} 
equipped  with  bhe  metric  (2.18).  Observe  that 

(2.20)  {u  }  c  X  and  lu  -  ul  .  +0  **>  u  e  X  and  d(u  ,u)  ♦  0 

n  T  n  l’cO)  T 

This  is  the  standard  remark  that  weakening  a  metric  of  a  compact  metric  space  produces  the 
same  topology. 

To  study  the  large  time  behaviour  of  u(t,uQ)  we  introduce  its  w-limit  set: 

y<uA)  •  {w  e  X  j  u(t  ,u„)  ♦  w  in  X  for  some  sequence  {t  }  with  t  +  •  as  n  +  «}. 

0  n  0  n  n 

We  collect  some  basic  remarks. 

Proposition  16.  Let  hypotheses  H  hold.  Then 

( i)  is  a  precompact  subset  of  X  for  x  >  0. 

(ii)  u(  *,u0)  e  C( (0,-)  :  X). 

(iii)  wtu^)  is  nonempty  and  connected  in  X. 

(iv)  If  we  uitUg),  then  u(t,w)  e  o)(uQ )  for  t  >  0. 


Proof .  Since  C  XT  w^ch  compact,  (i)  follows.  By  (2.20),  t  ♦  u(t,uQ)  is 

continuous  into  on  t  >  t  if  and  only  if  it  is  continuous  into  L^(f2),  whence  we 

have  (ii).  The  assertion  (iii)  follows  at  once  from  <i)  and  (ii).  For  (iv),  we  use  that 
u(t+tn'V  *  u(  t,u(tn,uQ ) )  so  if  u<tn,uQ)  ♦  w  in  X  (and  so  in  L^fl))  we  have 
u(t+tn,u0)  +  u(t,w)  in  L1  (and  hence  in  X)  by  Theorem  12.  Thus  u(t,w)  e  o)(uQ). 

Next  consider  the  function  V  s  X  -►  R  given  by 

V(c)  -  fa  (V2  (?V2  -  mF(?))dx 


where 


F(r)  *  pm  1  f  (  p)  dp  . 

Clearly  V  :  X  *  R  is  continuous * 

Lemma  1 7 .  Let  hypotheses  H  be  satisfied.  Then 

m+ 1 

(u  2  )t  e  s  L2(f)))  and 

m+1 

— JnUu  2  L  ^  +  V  ( u  ( t ,  u  ) )  <  V(u(s,u  ))  for  t  >  s  >  0  . 

(m+1  ) 

The  proof  is  postponed  briefly  while  we  establish  the  next  result. 

Theorem  18.  Let  hypotheses  H  be  satisfied.  Then  (o(uQ)  c  E. 

Proof .  By  Lemma  17  t  -►  V(u(t,uQ))  is  nonincreasing  on  t  >  0.  Since  V  is  continuous 
on  X  this  implies 


V(w) 


inf  V(u(t,uJ)  =  V 

_  U  00 

t>n 


for  w  e  (u(uj 
0 
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By  the  invariance  of  w(u  ),  V(u(t,w))  *  V  for  t  >  0,  w  e  a>( u _ )  -  Combining  this  with 

0  “  0 

m+1 

2 

Lemma  17  we  deduce  that  (u(t,w)  )  s  0  and  thus  u(t,w)  =  w.  The  definition  of  a 

solution  of  Problem  I  then  implies 


L(w%  +  f(w)  *)  «  0 
J  ft  xx 

whenever  V>e  C2(ft),  >  0  and  v?  (±L)  «  0*  But  this  implies  +  f(w)  *  0  in  P' 

Since  w  6  L°*  and  f  is  Lipschitz,  the  equation  holds  classically  (i.e.,  wm  e  C2(ft)) 
and  w  *  0  at  ±L.  Hence  w  e  E(L). 

OB 

Proof  of  Lemma  17.  We  go  back  into  the  existence  proof  where,  assuming  u^  e  CQ(ft),  we 

constructed  smooth  solutions  u  of  approximate  problems  I  •  Putting  t  -  0  in  (2.9) 

e  e 

and  letting  e  tend  to  zero  through  the  sequence  as  in  that  proof  yields 


t/o  fa  {(u  2  k}2  +  v<u(t'V’  <  v<u0> 


We  note  also  that  since  u™  -  em  e  ( Q)  for  t  >  0  and  (2.10)  holds,  we  have 

um  €  H^(ft)  for  t  >  0.  To  deduce  (2.21)  for  general  u„  e  H^(ft),  choose  (u„  }  c  c*(ft) 

0  00  On  0 

so  that  urt  +  urt  in  X  as  n  -►  ».  Writing  (2.21)  for  u  *  u(t,un  1  in  place  of  u 
On  0  n  un 

and  letting  n  >  •  establishes  the  inequality.  The  lemma  now  follows  if  we  show 
u(s,uQ)m  e  H^fft)  for  s  >  0.  But  we  know  u(s,u0)m  6  Hg(ft)  for  uQ  e  H^(ft)  and  that 
(u(s,u0)m)x  is  bounded  in  L2(ft)  (even  L°°(ft))  independently  of  uQ,  0  <  u  <  1,  whence 
the  result. 

Remarks  on  Part  II 

Section  3.  The  reader  should  notice  that  the  whole  development  of  this  section  is  valid  if 
ft  is  a  domain  in  rather  than  an  interval  of  R.  For  this  one  replaces  ri(u)^  by 


An(u)  (or  En(u)  where  E  is  a  suitable  elliptic  operator)  and  modifies  the  statement 


of  the  boundary  conditions  in  Problem  III.  The  estimate  (2.1)  remains  valid  as  stated. 
Various  of  the  ideas  used  in  this  proof  occur  for  example,  in  (14),  (15),  (17).  The  paper 
(19)  is  an  early  example  of  iJ-type  estimates. 

It  is  known  in  nonlinear  semigroup  theory  that  Problem  III  (with  ft  c_  rN)  has  a 
unique  solution  in  the  sense  of  this  section  if  ^  -  0  and  g  e  •  Moreover,  for 

this  n  need  only  be  continuous  (not  locally  Lipschitz  continuous)  and  the  estimates  (2.1) 
and  (2.6)  are  valid  for  these  solutions  (with  ^  *  0).  However,  the  proof  goes  by  showing 
uniqueness  (without  establishing  the  estimates)  and  obtaining  the  estimates  in  the 
construction  of  the  solutions.  See  [6]  and  its  references  concerning  the  uniqueness.  See, 
e.g.,  (1]f  (8),  [12]  concerning  the  semigroup  theory.  One  will  not  find  the  claims  above 

presented  clearly  in  these  sources,  they  are  (true)  folk-lore.  Moreover,  the  semigroup 

N  1 

theory  provides  solutions  to  Problem  III  (with  ft  c  X #  ^  *  0)  if  g  is  merely  L1  and 

Uq  is  merely  L1  and  n  need  not  be  a  function  but  a  graph.  One  can  also  take 

Uq  e  H  ^  ( ft) ,  f  e  (0 , T  :  H  1  ( C2 ) )  •  However,  one  does  not  then  use  the  above  notion  of 

solution  -  in  particular,  u  need  not  be  bounded.  See  [5]  concerning  the  H”1  theory. 
Similar  remarks  pertain  to  Problem  I,  although  now  it  is  the  L1  semigroup  theory  (rather 
than  H~^)  which  should  be  used. 

Section  4.  We  have  given  the  quickest  existence  proof  suitable  for  our  purpose.  It  is 
rather  standard  and  restrictive  in  that  it  requires  uQ  >  0.  Alternatives  are  provided  by 

the  semigroup  theory  (see  above),  but  this  is  clumsy  as  regards  approximations  by  smooth 

functions.  To  allow  Ug  to  change  sign  and  n  to  be  less  regular,  one  can  approximate 
n  by  smooth  end  regularize  the  equation  by  u^  -  AfrMu)  +  eu)  +  f(u).  (Again  the 

proofs  in  this  section  work  for  ft  c  x  )• 

The  argument  for  taking  the  limit  in  this  section  works  in  essence  if  luj1*"^  is 
replaced  by  n(u)  where  n  is  strictly  increasing  and  either  n  or  n  1  is  Lipschitz 
continuous.  For  general  n,  the  semigroup  theory  is  best.  With  it  one  has  that 
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-  An(u)  +  f(u) 


in  ftp  , 


1  3fl 


u|t-o  “  uo 


t  >  0  , 

in  fl  , 


has  a  solution  (in  the  semigroup  sense)  under  very  mild  conditions  which  we  do  not  detail 
here.  See  also  [11],  [3]. 


Section  5 .  The  relevant  reference  here  has  been  given  -  [6].  This  section  depends  on 
N  *  1  in  that  an  estimate  of  Aum  in  the  space  of  measures  provides  compactness  of  Tu™ 
in  only  if  N  *  1.  Moreover,  we  use  strongly  the  special  form  |u|m-1u  of  the 

nonlinearity.  The  only  related  work  on  more  general  nonlinearities  we  know  of  is  (9), 

[to]. 


Section  6«  Owing  to  the  remarks  concerning  Section  5,  these  arguments  do  not  adapt  to 
N  >  1.  Stabilization  results  in  higher  dimensions  remain  an  interesting  open  problem. 


-33- 


UMt' 


REFERENCES 


[1]  Barbu,  V.,  Nonlinear  Semigroups  and  Differential  Equations  in  Banach  Spaces, 

Noordhoff  International  Publishing  Co.,  Leyden  <1976). 

[2]  Bear,  J.,  Dynamics  of  flulda  In  porous  media,  American  Elzevier,  New  York,  1972. 

[3]  Benilan,  Ph.,  Sur  un  problems  d' evolution  non  monotone  dans  Publications 

Mathematiques  du  Besancon. 

[4]  Benilan,  Ph.  and  M.  G.  Crandall,  Regularizing  effects  of  homogeneous  evolution 
equations,  Mathematics  Research  Center  TSR  #2076,  University  of  Wisconsin-Madison, 
1980. 

[5]  Brezia,  H. ,  Monotonicity  methods  in  Hilbert  spaces  and  some  applications  to  nonlinear 
partial  differential  equations.  Contributions  to  Nonlinear  Functional  Analysis,  E. 
Zarantonello,  ed. ,  Academic  Press,  New  York  1971. 

[6]  Brezis,  H.  and  M.  G.  Crandall,  Uniqueness  of  solutions  of  the  in it ial- value  problem 
for  u^  -  A^(u)  *  0,  J.  Math.  Pures  et  Appl.  58  (1979)  153-163. 

[7]  Chafee,  N.  and  E.  Infante,  A  bifurcation  problem  for  a  nonlinear  parabolic  equation, 
Applicable  Analysis  4  (1974),  17-37. 

{8]  Crandall,  M.  G. ,  An  introduction  to  evolution  governed  by  accretive  operators, 

Dynamical  Systems  -  An  International  Symposium,  L.  Cesari,  J.  Hale,  J.  LaSalle,  eds., 
Academic  Press,  New  York,  1976,  pp.  131-165. 

[91  Crandall,  M.  G.  and  M.  Pierre,  Regularizing  effects  for  u^  «  A^(u) ,  Mathematics 
Research  Center  TSR  #2166,  University  of  Wisconsin-Madison,  1981. 

[101  Crandall,  M.  G.  and  M.  Pierre,  Regularizing  effects  for  -  A<^(  u)  in  L1, 
Mathematics  Research  Center  TSR  #2187,  University  of  Wisconsin-Madison,  1981. 

[11]  Evans,  L.  C.,  Differentiability  of  a  nonlinear  semigroup  in  L1,  J.  Math.  Anal.  & 
Appli.,  60  (1977),  703-715. 

[12]  Evans,  L.  C.,  Application  of  nonlinear  semigroup  theory  to  certain  partial 
differential  equations.  Nonlinear  Evolution  Equations,  M.  G.  Crandall,  ed.,  Academic 
Press,  New  York  1978. 

[13]  Gurtin,  M.  E.  and  R.  C.  MacCamy,  On  the  diffusion  of  biological  populations, 
Mathematical  Biosciences  33  (1977),  35-49. 

[14]  Kamin,  S.,  Source-type  solutions  for  equations  of  nonstationary  filtration,  J.  Math. 
Anal.  Appl.  63  (1978). 

[15]  Ladyzenska ja,  O.  A.,  V.  A.  Solonnikov  and  N.  N.  Ural'ceva,  Linear  and  Qiasilinear 
Equations  of  Parabolic  Type,  Transl.  of  Math.  Monographs  Vol.  23,  American 
Mathematical  Society,  Providence,  1968. 

[16]  Okubo,  A.,  Diffusion  and  Ecological  Problems:  Mathematical  Models,  Biomathematics, 
Vol.  10,  Springer-Verlag,  Berl in- Heidelberg- New  York,  1980. 


[17]  Pierre#  M. ,  Uniqueness  of  the  solutions  of  ufc  -  <V(u)  *  0  with  initial  datum  a 

measure#  Mathematics  Research  Center  TSR  #2171,  University  of  Wisconain-Madison , 
1981. 

[18]  Smoller#  J.  and  J.  Wasserman#  Global  bifurcations  steady-state  solutions,  J.  Diff. 
Bqns.  39  (1981),  269-290. 

[19]  Vol'pert,  A.  I.  and  S.  I.  Hudjaev,  Cauchy's  problem  for  degenerate  second  order 
quasi linear  parabolic  equations.  Math.  USSR  Sbornik,  7  (1969),  365-387. 


DA/MGC/LAP/jvs 


SECURITY  CLASSIFICATION  OF  This  PAGE  (When  Data  Entered ) 


REPORT  DOCUMENTATION  PAGE 


1.  REPORT  NUMBER 

#2220  v 


4.  TITLE  (end  Subtitle) 

Stabilization  of  Solutions  of  a  Degenerate 
Nonlinear  Diffusion  Problem 


7.  AUTHORf*) 

Donald 'Aronson ,  Michael  G. 1  Crandall  and 
L.  A.  ,-Peletier 


#.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Mathematics  Research  Center,  University  of 
610  Walnut  Street  Wisconsin 

Madison.  Wisconsin  53706 


11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

U.  S.  Army  Research  Office 
P.O.  Box  12211 

Research  Triangle  Park,  North  Carolina  27709 


.  MONITORING  AGENCY  NAME  A  ADDRESS (it  different  from  Controlling  Office) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


5.  TYPE  OF  REPORT”®  PERIOD  COVERED 

Summary  Report  -  no  specific 
reporting  period 


6.  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUMBER^*) 

'  DAAG29-80-C-0041* 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 

Work  Unit  Number  1  - 
Applied  Analysis 


12.  REPORT  DATE 

May  1981 


NUMBER  OF  PAGES 

35 


IS.  SECURITY  CLASS,  (of  thia  report) 

UNCLASSIFIED 


t5«.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  thia  Report) 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abatrect  entered  In  Block  20,  if  different  from  Report) 


19.  KEY  WOROS  (Continue  on  reeeree  aide  if  nmcaaaary  and  identify  by  block  number) 


nonlinear  diffusion,  biological  models,  asymptotic  behaviour,  equilibrium 
solutions,  flow  in  porous  media 


20.  ABSTRACT  (Continue  on  reeeram  aide  if  necaaasry  and  Identity  by  block  number) 

The  central  results  concern  the  initial* value  problem  for  u  =  (u™)  + 

t  xx 

u(l-u)  (u-a)  ,  -L  ±  x  <  L  and  t  j>  0,  under  Dirichlet  conditions  at  x  -  ±L. 

Here  m  >  1  and  0  <  a  <  (m+1) /(m+3)  -  The  equilibrium  solutions  of  this  problem 
are  determined  for  each  L  >  0  and  it  is  shown  that  the  u>-limit  set  gj(u0)  of 

an  initial  datum  with  values  in  10,1)  consists  of  a  connected  set  of  equi¬ 

libria.  This  is  used  to  determine  some  domains  of  attraction  of  isolated 


1  jAN^S  1473  COITION  OF  1  NOV  6*  IS  OBSOLETE 


UNCLASSIFIED 

ftCCUfttTY  CLASSIFICATION  OF  THIS  PAGE  (When  t>ef  Entered) 


ABSTRACT  (continued) 


equilibria.  A  novel  feature  of  the  results  is  that  for  large  L  there  are 
multiple  parameter  families  of  equilibria. 

A  second  part  of  the  paper  gives  a  self-contained  development  of  existence, 
uniaueness,  maximum  principles,  and  continuous  dependence  on  data  for  more 
general  equations  u^  =  n(u)^  +  f(u).  The  results  are  employed  in  proofs  of 

some  of  the  theorems  referred  to  above. 

Interest  in  these  questions  is  stimulated  by  the  occurrence  of  such  models 
in  science,  e.g.  in  fluid  flow  in  porous  media  and  biology* 


\ 


